Comment on "First-principles calculation of the superconducting transition in MgB2 

within the anisotropic Eliashberg formalism" 
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Choi et al. [Phys. Rev. B 66, 020513 (2002)] recently presented first principles calculations 
of the electron-phonon coupling and superconductivity in MgB2, emphasizing the importance of 
anisotropy and anharmonicity. We point out that (1) variation of the superconducting gap inside 
the a- or the 7r-bands can hardly be observed in real samples, and (2) taking the anisotropy of the 
Coulomb repulsion into account influences the size of the small gap, A-n-. 
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In a recent paper^i as well as in a follow-up paper,'^ 
Choi et al. presented an ah initio calculation of the su- 
perconducting transition and superconducting properties 
of MgB2. The important improvement over existing cal- 
culations was that they allowed the order parameter to 
vary freely over the Fermi surface, i.e. A = A (k), and at 
the same time took the anharmonicity into account. As 
a consequence, they had to compute the fully anisotropic 
electron-phonon interaction, A (k, k') , and to solve the 
corresponding Eliashberg equation. The Coulomb pseu- 
dopotcntial, /i* (k, k') , was assumed not to depend on 
k and k', and was treated as an adjustable parame- 
ter. First-generation ah initio calculations of the super- 
conducting transition and superconducting properties of 
MgB2 had assumed A to be constant and had therefore 
solved merely the isotropic Eliashberg equation.'^ More- 
over, anharmonicity had been neglected. It was soon 
pointed ouf* that the calculated electron-phonon cou- 
pling suggests that the gap on the two tt sheets of the 
Fermi surface is smaller than that on the two a sheets, 
and that anharmonicity is important. This led to the so- 
called two-band model. Ab initio calculations of the sec- 
ond generatio niiai^ allowed for two, and sometimes four 
gaps, A„, and thus had to compute A„„' , to estimate 
the anisotropy of ^*„/, and to solve the corresponding 
Eliashberg equations. 

Here we shall comment on (1) whether consequences 
of anisotropy beyond that of the two-band model may 
be observed and (2) whether at this level of detail Choi 
et a/.'s assumption of a uniform Coulomb repulsion is 
warranted. 

1. RefQ implies that there is a distribution of gaps 
within the a and the 7r-sheets, not only in the calcula- 
tions for perfectly clean MgB2, but also in the actual 
material; in other words, that the distribution of gaps 
shown in Fig. 2 of Ref 2 is observable. However, in the 
theory of anisotropic superconductivity it is known that 
any intraband nonuniformity of the order parameter is 
suppressed by strong intraband impurity scattering. It is 
not immediately obvious, though, when scattering should 
be considered strong in this connection. Since excitation 
gaps are not equal to the order parameters any more, one 



needs to compare individual densities of states (DOS), 
N{E), for the two a-bands (or the two 7r-bands) with 
each other for a given scattering strength, and check 
whether \Nai—Na2\ ^ \Nai+Na2\- The relevant expres- 
sion can be found in Ref 7. In the limit of large scattering 
rates, 7, one can derive an analytical expression for this 
criterion, namely 7 > yj (A) ■ i5A (Ref@), where (A) is 
the average order parameter, and JA is the variation of 
the order parameter over the Fermi surface in question. 
With the data from RefsQ0 for (A) and (5A, this gives 
characteristic scattering rates of respectively 2 and 1.5 
meV for the a and 7r-bands. Therefore, to observe 4 dis- 
tinct gaps in MgB2 one needs samples with scattering 
rates smaller than 2 meV, that is, with mean free paths 
beyond 1500 A. To observe gap variations beyond the 4- 
band model, far cleaner samples are needed. This is the 
reason why at most two distinct gaps have been observed 
in experiments. It is even surprising that the difference 
of 5 meV between the gaps of the a and the 7r-bands is 
not smeared out. This seems to be due to the inability 
of common impurities to couple between the disparate a 
and TT-band wavefunctionsp so that ^^tt ^ laa ^ 



2. For the Coulomb pseudopotcntial, Choi et al. used 
^*(k, k') ~ h*{llIc) = 0.12 (with the cut-off frequency 
ujc « bujp]^^) and stated that the superconducting prop- 
erties of MgB2 were not very sensitive to the choice 
of /i*((jJc)- This at first seems plausible, because the 
Coulomb pseudopotcntial enters the Eliashberg equation 
only in the combination A(k, k', 1^ — ly') — /i*(k, k'), and 
the A-distribution varies on the scale of ~1.8, ~0.3, and 
^0.2 for CTcr-, TTTT-, and crTr-scattering respectively [see Fig. 
3 of Ref Ij. Therefore, at most the tT7r-scattering can be 
influenced by anisotropy of /i*. We shall argue that the 
CTTT-interband Coulomb matrix elements are considerably 
smaller than the intraband matrix elements due to the 
very small overlap of the a- and 7r-band charge densities^ 
and that this is sufficient to infiuence the superconduct- 
ing properties, in particular the size of the small gap, 
A^. 

Choi et al. do not give the band-integrated values of 
their coupling constants, but by integrating Fig. 3 of 
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RefQ with the DOS-ratio N^/N„ = 1.37 according to: 

A„„,(0) ^ l^fellA(k,k',0)<5(£„-k') (1) 
k,k' " 

^ 5Z A(k, k',0)W^„-k' A^„' 



N 

k,k' 

for the phonon-mediated coupling of an electron in band 
n to all electrons in band n', we can map Choi et 
al.^s fully anisotropic model onto a two-gap model with 
Xaa = 0.78, Xan = 0.15, Ka = 0.11, and A^^ = 0.21. 
These A-values yield the mass-renormalization parame- 
ters in Fig. 2 of Ref.l!: m*/m-l = \„ = \„^+Xaa ~ 0.94 
and Att — Attct -I- A^ttt ~ 0.32. The total isotropic (thermo- 
dynamic) A = (iVcrAcr +iV7rA7r)/iV = 0.61, which of course 
is the same as the one given by Choi et al. Here, and 
in Eq. ^ is the DOS summed over all bands. With 
this two-gap model we have performed strong-coupling 
Eliashberg calculations in order to compare the results 
for Tc and the gaps with those resulting from the fully 
anisotropic treatment. For all four spectral functions we 
used the isotropic a'^F{uj) from Fig. 1 of RefQ scaled to 
produce the A-matrix given above. The /i*(a;c)-matrix 
is obtained from Eq. ^ with A(k, k',0) substituted by 
Choi et al.'s iJ.*{u!c)- The resulting Tc and the gaps are 
shown by dashed Hnes in Fig. 1 as functions of ^*{uJc)- 
At fi*{ujc) = 0.12, as used by Choi et al, we get Tc = 43 
K, Ao- = 7.2 meV, and A^r = 1.3 meV. The correspond- 
ing values quoted by Choi et al. are 39 K, 6.8 meV and 
1.8 meV. These differences are hardly due to intraband 
anisotropy, first of all because it can only increase Tc- 
Secondly, increasing the number of gaps from two to four 
in the Eliashberg equations, which should account for 
most of the anisotropy beyond the two gap model, we 
found rather small changes »iS 

If, on the other extreme, we assume that there is no 
Coulomb repulsion between the a and 7r-electrons, then 
the corresponding two-gap treatment gives the full lines 
in Fig. 1 and, hence, Tc — 38 K, Ao- = 6.5 meV, and 
Att — 1.8 meV for ii*{u)c) = 0.12, incidentally, rather 
close to the values quoted in Refs[l||3- If the magnitude 
of /i* in both calculations shown in Fig. 1 is adjusted to 
produce the same Tc of 39 K, the value of the lower gap 
changes from « 2 eV (diagonal) to w 0.4 eV (uniform). 

That uniform and diagonal Coulomb pseudopotentials 
yield different results is not surprising: The same to- 
tal Eliashberg /x* in the uniform case is distributed over 
intra- and interband terms so that the crCT-part of the 
pairing interaction suffers less than in the case of a di- 
agonal /I*. Xaa is more important for the critical tem- 
perature, and Acttt for generating A^r. For uniform /i*, 
therefore, the Tc and A^ are larger, and A^ is much 
smaller. 

Having demonstrated that the assumed structure 
of ^* matters for the details of the superconducting 
properties of MgB2ri the size of A^ in particular, 
let us finally estimate this structure from first prin- 
ciples. The unrenormalized /i is the matrix element 
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FIG. 1: Critical temperature and the values of the a and tt 
gaps at 1 K as function of the renormalized Coulomb pseu- 
dopotential, /i*(ajc), in two models: The uniform model where 
all matrix elements of the Coulomb repulsion are equal and 
the diagonal model where the interband matrix elements are 
zero. In both cases the normalization is chosen so as to pro- 
duce given values of ^*{uc) after proper summation over all 
bands. The two gaps obtained from /x*(tJc)'s giving T^s of 39 
K in the two models are connected vertically. 



(nk t,n — k llVcln'k' f , n' — k' J,) for scattering a Cooper 
pair from state |n'k') to state |nk) via a phonon with 
wave- vector k — k'. Inserting this matrix element in Eq. 
n instead of A(k, k',0) yields ^„„'. Here Vc{y,t') is the 
screened Coulomb interaction between the electrons, and 
since it has short range in good metals, it makes sense to 
take it proportional to the delta- function 5{v — r'). This 
leads to the following estimate: 

fi cx N [ |V^„k(r)|'|V^„'k'(r)|'rf'r 



/inn' cx iV„, / |V'(r)|^^ \^,{v)\l, (2) 



where |-0(r)|„ = I]k l'^nk(r)| 5{en'k)/Nn is the shape, 
normalized to 1 in the cell or the crystal, of the electron 
density of band n at the Fermi level. These cr and tt 
densities are shown in Fig. 2, and they yield for the ratios 
of the integrals in Eq[21 



3.0 : 1.8 : 1. 



(3) 



These ratios reflect the facts that the cr-density is more 
compact than the 7r-density, and that the overlap of these 
two densities is small. Note that the exceptional small- 
ncss of the interband impurity scattering in MgB2 is due 
not only to this difference in charge density, but also to 
a disparity of the a and tt wave functions. 

From Eqs. |21lllwe get: fj,cra-fJ'TnT-fJ-a-w-f^7ra—'^-^-'2.Q.l.4::l. 
Now, any anisotropy in the bare pseudopotential is fur- 
ther enhanced in the renormalized fi*: In the one-band 
case fi is renormalized as ^*{uJc) — ^j./ [I + ^IniW/uc)], 
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FIG. 2: Constant-density contour for the normalized a (or- 
ange) and TT (green) electron-densities, ('(/'(r)]^ and |?/)(r)|^ , 
at the Fermi level. 



where W is a characteristic electronic energy of the order 
of the bandwidth or plasma frequency. For the multi- 
band case, this is a matrix equation with W being a 
diagonal matrix with elements Wn- Assuming for sim- 
plicity that fiaa — fJ-TTTT — A^cy-K with A > 1, and that 
^crcr log(VFcr/wc) = log( /tjc ) = L, ouc obtalus: 
A* ^ A + [A - A-^)L. For MgBa, L - 0.5 - 1 and 
A 2.3, so that A* ~ 3 — 4, which is very different from 
the uniform /i. 



In conclusion: Any difference between the results of 
the fully anisotropic Eliashberg formalism and those of 
the two-gap formalism will hardly be observable in real 
MgB2-samples. On the other hand, the anisotropy of the 
Coulomb pseudopotential is likely to have an observable 
effect on the size of the small gap, A^r. 



The authors thank W.E. Pickett for numerous helpful 
discussions, and for critical reading of the manuscript. 
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We have verified this analytical criterion numerically using 
some simplified models for the electron-phonon coupling 
matrix. 
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° This is not surprising, given that the variation of the 
calculated order parameter;- as well as of the mass 
renormalization,^ within each band is at most 10%. 

^ The problem of proper structure of the /i*-matrix is specific 
for multiband superconductivity with large gap disparity. 
When the pairing interaction is fully uniform, the n* ma- 
trix can have any structure as long as it is normalized to 
the total fi* . This has no effect on either Tc or gap-ratios. 



